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Abstract. I describe a class of oscillating bounce solutions to the Euclidean field equations for gravity coupled to a scalar 
field theory with multiple vacua. I discuss their implications for vacuum tunneling transitions and for elucidating the thermal 
nature of de Sitter spacetime. 


INTRODUCTION 


The problem of vacuum tunneling in de Sitter spacetime 
has recently acquired renewed relevance. In part, this is 
due to developments in string theory, which suggest that 
vacuum tunneling may be of relevance for understand¬ 
ing transitions between the various potential vacua that 
populate the string theory landscape. But, it is also of in¬ 
terest for the light that it can shed on the nature of de 
Sitter ^acetime. In this talk I will describe some recent 
work llll] with Jim Hackworth in which we explored some 
aspects of the subject that have received relatively little 
attention. 

De Sitter spacetime is the solution to Einstein’s equa¬ 
tions when there is a constant positive vacuum energy 
density Vyac, but no other source. Globally, it can be rep- 
resented as the hyperboloidx +y +z +w —v = 
in a flat five-dimensional space with metric ds^ = dx^ + 
dy^ + dz^ + dw^ — dv^, where 
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The surfaces of constant v are three-spheres, with the 
sphere of minimum radius, occurring at v = 0. 

However, the special role played by this surface is 
illusory. De Sitter spacetime is homogeneous, and a 
spacelike three-sphere of minimum radius can be drawn 
through any point. 

An important property of de Sitter spacetime is the 
existence of horizons. Just as for the case of a black hole, 
the existence of a horizon gives rise to thermal radiation, 
characterized by a temperature 


Tis=H/2n. (2) 


However, there are important differences from the black 
hole case. A black hole horizon has a definite loca¬ 
tion, independent of the observer. Further, although an 


observer’s motion affects how the thermal radiation is 
perceived, the radiation has an unambiguous, observer- 
independent consequence — after a finite time, the black 
hole evaporates. By contrast, the location of the de Sit¬ 
ter horizon varies from observer to observer. Although 
comoving observers detect thermal radiation with a tem¬ 
perature Tjjs, this radiation does not in any sense cause 
the de Sitter spacetime to evaporate. As we will see, tun¬ 
neling between different de Sitter vacua provides further 
insight into the thermal nature of de Sitter spacetime and 
the meaning of T^s- 

Of course, the relevance of de Sitter spacetime to our 
Universe comes from the fact that in the far past, during 
the inflationary era, and in the far future, if the dark 
energy truly corresponds to a cosmological constant, the 
Universe approximates a portion of de Sitter spacetime. 
The underlying assumption is that results derived in the 
context of the full de Sitter spacetime are applicable to 
a region that is approximately de Sitter over a spacetime 
volume large compared to 
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FIGURE 1. The potential for a typical theory with a false 
vacuum 




VACUUM DECAY IN FUAT SPACETIME 


The prototypical example for studying vacuum decay 
is a scalar field theory with a potential, such as that 
shown in Fig. ^ that has both an absolute minimum 
(the “true vacuum”) and a higher local minimum (the 
“false vacuum”). For the purposes of this talk I will 
assume that y(^) > 0, so that both vacua correspond to 
de Sitter spacetimes. The false vacuum is a metastable 
state that decays by quantum mechanical tunneling. It 
must be kept in mind, however, that the tunneling is not 
from a homogeneous false vacuum to a homogeneous 
true vacuum, as might be suggested by the plot of V (0). 
Rather, the decay proceeds via bubble nucleation, with 
tunneling being from the homogeneous false vacuum to 
a configuration containing a bubble of (approximate) true 
vacuum embedded in a false vacuum background. After 
nucleation, the bubble expands, a classically allowed 
process. 

I will begin the discussion of this process by recalling 
the simplest case of quantum tunneling, that of a point 
particle of mass m tunneling through a one-dimensional 
potential energy barrier U (q) from a initial point to 
a point q'fln on the other side of the barrier. The WKB 
approximation gives a tunneling rate proportional to 
where 

B = 2 f dqy^2m[U(q) — E\. (3) 

This result can be generalized to the case of a multi¬ 
dimensional system with coordinates qi,q 2 ,...,qN- 
Given an initial point one considers paths qj{s) 
that start at qj"^ and end at some point q^" on the 
opposite side of the barrier. Each such path defines a 
one-dimensional tunneling integral B. The WKB tunnel¬ 
ing exponent is obtained from the path that minimizes 
this integral As a bonus, this minimization process 
also determines the optimal exit point from the barrier. 

By manipulations analogous to those used in classi¬ 
cal mechanics (but with some signs changed), this mini¬ 
mization problem can be recast as the problem of finding 
a stationary point of the Euclidean action 
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One is thus led to solve the Euclidean equations of mo¬ 
tion 


0 = m 




du 


(5) 


The boundary conditions are that qj{Ti„if) = g'™* and 
(because the kinetic energy vanishes at the point where 
the particle emerges from the barrier) that dqj/dx = 0 
at Tfin- The vanishing of dqj/dx at the endpoint implies 
that the solution can be extended back, in a “r-reversed” 


fashion, to give a solution that runs from to q^^ and 
back again to This solution is known as a “bounce”, 
and the tunneling exponent is given by 


B = 





= (bounce) — (false vacuum), 


( 6 ) 


with the factor of 2 in Eq. being absorbed by the 
doubling of the path. It is essential to remember that x is 
not in any sense a time, but merely one of many possible 
parameterizations of the optimal tunnelin^ath. 

The translation of this to field theory is straight¬ 
forward; The coordinates qj become the field variables 
())(x), and the path qj{x) becomes a series of three- 
dimensional field configurations ())(x,t). The Euclidean 
action is 



dxd^x 
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and so one must solve 


d^ij) 

dx'^ 




dV 
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The boundary conditions are that the path must start 
at the homogeneous false vacuum configuration, with 
())(x,Tinit) = 0fv, and that d(j)/dx = 0 at Tfln for all x. 
(Because is a minimum of the potential, it turns out 
that Tinit = —°°.) A three-dimensional slice through the 
solution at Tfin gives the most likely field configuration 
for the nucleated true vacuum bubble. This configura¬ 
tion, ^(x,Tfin), gives the initial condition for the sub¬ 
sequent real-time evolution of the bubble. As with the 
single-particle case, a T-reflected solution is convention¬ 
ally added to give a full bounce. 

Despite the fact that the spatial coordinates x and the 
path parameter x have very different physical meanings, 
there is a remarkable mathematical symmetry in how 
they enter. This suggests looking for solutions that have 
an SO(4) symmetry; i.e., solutions for which 0 is a 
function of only s = vx^ + z^. Eor such solutions, the 
field equation reduces to 

3 d^ dV 

ds^ s ds dcj) 

The boundary conditions are 

= 0 , = ( 10 ) 

i=0 

where the first follows from the requirement that the 
solution be nonsingular at the origin, and the second 


ds 












ensures both that a spatial slice at T = corresponds to 
the initial state, and that the slices at finite T have finite 
energy relative to the initial state. Note that while ^(0) is 
required to be on the true vacuum side of the barrier, it is 
not equal to (although it may be close to) 0tv 

Although the tunneling exponent is readily obtained 
from the WKB approach, the prefactor, including (in 
principle) higher order corrections, is most easily calcu¬ 
lated from a path integral approach 01 ■ The basic idea 
is to view the false vacuum as a metastable state with 
a complex energy, with the imaginary part of the energy 
density yielding the decay rate per unit volume. The false 
vacuum energy is obtained by noting that for large ^ 

I{3r) = j^ , ( 11 ) 

where the path integral is over configurations with 

^{x,T = ±£r/2) = ^fy. 

This path integral can be calculated by summing the 
contributions from the various stationary points, each of 
which gives a factor of (det Here S" is the 

functional second derivative of the action, evaluated at 
the stationary point; i.e., the product of the frequencies 
of the normal modes. The first stationary point, a homo¬ 
geneous false vacuum configuration with 0(x,t) = 0fv 
everywhere, gives a contribution where the real 

prefactor A includes the (properly renormalized) deter¬ 
minant factor and 5fv = V{(j)fv) Here Y denotes the 
volume of space and is understood to be taken to infinity 
at the end of the calculation. 

The next stationary point is the bounce solution to 
Eq. The calculation of the determinant factor here is 
complicated by the fact that 5"(^ounce) has one negative 
and four zero eigenvalues. The former implies a factor of 
i, which 1 will display explicitly. The latter require the 
introduction of collective coordinates; integrating over 
these gives a factor of corresponding to the fact 
that the bounce can be centered anywhere in the four¬ 
dimensional Euclidean space. 

Einally, the approximate stationary points corre¬ 
sponding to multibounce solutions also contribute, 
with the n-bounce contribution including a factor of 
{i7'yY/n\ from integrating over the positions of n 
identical bounces. Putting all this together gives a result 
that can be written as 


I[3r) = Ae ^^''+iy^Je -|- 
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point is that it is real. Extracting the energy density from 
the exponent in Eq. o gives an imaginary part that is 
proportional to fT, corresponding to the fact that a bubble 
can nucleate anywhere. The quantity we actually want is 
the nucleation rate per unit volume. 
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The path integral approach provides the vehicle for 
extending the calculation to finite temperature T, 
with the path integral over configurations extending from 
T = —ooto T = oo replaced by one over over configu¬ 
rations that are periodic in T with periodicity I/T. At 
low temperature, where l/T is larger than the character¬ 
istic radius of the four-dimensional bounce, there is little 
change from the zero-temperature nucleation rate. How¬ 
ever, in the high-temperature regime where l/T is much 
smaller than this characteristic radius the path integral 
is dominated by configurations that are constant in T. A 
spatial slice at fixed T gives a configuration, with total 
energy Ecrit, that contains a single critical bubble. The 
exponent in the nucleation rate takes the thermal form 



(14) 


Note that, in contrast to the zero temperature case, there 
is no spatial slice corresponding to the initial state. Only 
through the boundary conditions at spatial infinity does 
the bounce solution give an indication of the initial con¬ 
ditions. 


ADDING GRAVITY 

Coleman and De Luccia II argued that the effects of 
gravity on vacuum decay could be obtained by adding an 
Einstein-Hilbert term to the Euclidean action and then 
seeking bounce solutions of the resulting field equations; 
as before, the tunneling exponent would be obtained 
from the difference between the actions of the bounce 
and the homogeneous initial state. Their treatment did 
not include the calculation of the prefactor, an issue that 
remains poorly understood. 

If one assumes 0(4) symmetry, as in the flat spacetime 
case, the metric can be written as 

ds^ = + (15) 

where dO.^ is the metric on the three-sphere, and the 
scalar field depends only on The Euclidean action 
becomes 


Here J includes both determinant and Jacobean factors, 
with the latter arising from the introduction of the col¬ 
lective coordinates; for present purposes, the important 


Se = 27f 




U^ + v 






3A/pi / 2 ■■ -2 \ 

-^{p p+pp -p) 


, (16) 


with dots denoting derivatives with respect to The 
Euclidean field equations are 


and 
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(18) 


One can show that if E((/)) is everywhere positive, as 
1 am assuming here, then p(^) has two zeros and the 
Euclidean space is topologically a four-sphere. One of 
the zeros of p can be chosen to lie at ^ = 0, while 
the other is located at some value ^max- Requiring the 
scalar field to be nonsingular then imposes the boundary 
conditions 


dtp 

d^ 



d(j) 

d^ 
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(19) 


The symmetry of these boundary conditions should be 
contrasted with the flat space boundary conditions of 
Eq. (ITO). Note that there is no requirement that scalar 
field ever achieve either of its vacuum values, although 
|0(^max) — 0fv| is typically exponentially small in cases 
where gravitational effects are small. 

Somewhat surprisingly, the Euclidean solution corre¬ 
sponding to a homogeneous false vacuum is not an infi¬ 
nite space, but rather a four-sphere of radius 


^fv = 


S7tV{(pfy) ■ 


Its Euclidean action is 
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( 21 ) 


If the parameters of the theory are such that the charac¬ 
teristic radius of the flat space bounce is much less than 
then the curved space bounce will be roughly as 
illustrated in Eig. with the small region near ^ — 0 
corresponding to the true vacuum region of the flat space 
bounce, the equatorial slice giving the optimal configu¬ 
ration for emerging from the potential barrier, and a slice 
such as that indicated by the lower dotted line roughly 
corresponding to the state of the system before the tun¬ 
neling process. A bounce solution such as this yields a 
nucleation rate that only differs only slightly from the 
flat space result. 

On the other hand, there are choices of parameters that 
give a bounce solution similar to that indicated in Eig. Eh, 
with a true vacuum region that occupies a significant 


fraction of the Euclidean space. In this case, there is no 
slice that even roughly approximates the initial state, sug¬ 
gesting that one should view this as more analogous to a 
thermal transition in flat space than to zero-temperature 
quantum mechanical tunneling. 

Indeed, for a bounce such as this the true and false 
vacuum regions can perhaps be viewed as being on a 
similar footing, so that the bounce can describe either 
the nucleation of a true vacuum bubble in a region of 
false vacuum, or the nucleation of a false vacuum bubble 
in a true vacuum region @1. The rate for the former case 
would be 


Efv^tv ^ exp{-[5E(bounce) -5E(fv)]} , (22) 

while for the latter. 


Rtv^fv ^ exp{-[5E(bounce) -5E(tv)]} . (23) 


The ratio of these is 


Rtv^fv 

rfv^tv 


exp{5E(tv) -5E(fv)} 
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(24) 


If y(0fv) — y(0tv) y(0fv)j the geometry of space is 
roughly the same in the two vacua, and we can sensibly 
ask about the relative volumes of space occupied by the 
false and true vacua. In the steady state, this will be 
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.(25) 


The last line of this equation, which gives the ratio as 
the exponential of an energy difference divided by the 
de Sitter temperature, is quite suggestive of a thermal 
interpretation of tunneling in this regime. 

It is not hard to show that the flat space Euclidean 
field equations always have a bounce solution. This is no 
longer true when gravity is included, as we will see more 
explicitly below. However, Eqs. (inii and Cl always 
have a homogeneous Hawking-Moss @] solution that is 
that is qualitatively quite different from the flat space 
bounce. Here (j) is identically equal to its value 0top at the 
top of the barrier, while Euclidean space is a four-sphere 

of radius 3Mpj/87ry(0top). Erom this solution 

one infers a nucleation rate 


f 3 3 ) 

1 8y((^)fv)^8y(0,op)j ■ 


( 26 ) 


from the false vacuum 


















Final state 


FIGURE 2. Schematic illustration of a Coleman-De Luccia bounce solution in two limiting regimes. In both, (j) is near its true 
vacuum value in the region to the right of the dashed arc, while on the left side it is near the false vacuum. In both cases the 
equatorial slice denotes a three-sphere corresponding to the spatial hypersurface on which the bubble nucleates. The lower dashed 
line in (a) represents a three-sphere indicative of the initial false vacuum state; this has no analogue in the regime illustrated in (b). 


OTHER TYPES OF BOUNCES? 


Given the existence of the Hawking-Moss solution, it 
is natural to inquire whether the inclusion of gravity 
allows any other new classes of Euclidean solutions. In 
particular, might there be “oscillating bounce” solutions 
in which (j) crosses the potential barrier not once, but 
rather k> I times, between ^ = 0 and ^ = £max? 

There can indeed be such solutions lUOl] . In exam¬ 
ining their properties, we focussed on the case where 
l^(0top) — T(0tv) ‘C y(0tv)- This simplifies the calcula¬ 
tions considerably, but does not seem to be essential for 
our final conclusions. With this assumption, the metric 
is, to a first approximation, that of a four-sphere of fixed 
radius and = 7t/H. We then only need to solve 
the scalar field Eq. (O- Defining y — , we can write 

this as 


d^<b d<b I dV 

—jT + 3 coty— = —=r — . 
dy^ dy d(j) 


(27) 


It is convenient to start by first examining “small am¬ 
plitude” solutions in which 0 (0) and ^ (tt) are both close 
to ^op- Let us assume that near the top of the barrier V 
can be expanded as' 


V(<?»)=V(<?>top) 


—^(0 —0top)^H——^top)"^H- 

(28) 


^ The omission of cubic terms here is only to simplify the algebra. 
There is no difficulty, and little qualitative change, in including such 
terms. The details are given in Ref. Q|]. 


with 




|V"(0top)| 

i/2 


(29) 


Keeping only terms linear in (0 — 0top) in Eq. (I17> gives 
d^(j) d(j) 


whose general solution is 


^(y) - 0top (cosy) (cosy), (31) 

where Cd and DJ, are Gegenbauer functions of the 
first and second kind and a(a -f 3) = j3. The vanishing 
of dcj)/dy at y = 0 implies that B = 0; the analogous 
condition at y = ;r is satisfied only if a is an integer, in 
which case Cd is a polynomial. 

While the linearized equation only has solutions for 
special values of j3, this condition is relaxed when the 
nonlinear terms are included. Eurthermore, the nonlinear 
terms fix the amplitude of the oscillations, which is com¬ 
pletely undetermined at the linear level. The problem can 
be analyzed by an approach similar to that used to treat 
the anharmonic oscillator. Any function with dcjt/dy van¬ 
ishing at both y = 0 and y = 7t can be expanded as 


^ (y) = + (y) • ( 32 ) 

Vl'^l M=0 

Substituting this into Eq. o and keeping terms up to 
cubic order in (0 — (jtiop) gives 


0 = E cli\cosy)\[P-M{M + 3)]AM 


M=0 













FIGURE 3. 


Bounce solutions for a scalar potential Q|] with cubic and quartic interactions and p = 70.03. 
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where the qijK\M arise from expanding products of three 
Gegenbauer polynomials. Requiring that the quantities 
multiplying each of the separately vanish yields an 
infinite set of coupled equations. These simplify, how¬ 
ever, if |A| = 1)3 — A^(A^ + 3)| <C 1 for some N. In this 
case, one coefficient. An, is much greater than all the oth¬ 
ers. The M — N term in Eq. 03 then gives (to leading 
order) 

An = ±\ -TjT-, (34) 

y sgn{A)qNNN-N 

where qNNN;N > 0- 

If A > 0, Eq. (I34t only gives a real value of An if )3 > 
N{N -I- 3). As j3 is increased through this critical value, 
two solutions appear. These are essentially small oscil¬ 
lations about the Hawking-Moss solution, with (j) (0) ~ 

(lkop±ANCl^^il) and (PiTt) « (lkop±ANCf/\-l). Be¬ 
tween these endpoints, (j) crosses the top of the barrier 


N times. If N is even, the two solutions are physically 
distinct, with one having (j) on the true vacuum side of 
the barrier at both endpoints, and the other having both 
endpoint values on the false vacuum side. If N is odd, the 
two solutions are just “y-reversed” images of each other. 

As j3 is increased further, the endpoints move down 
the sides of the barrier, until eventually the small am¬ 
plitude approximation breaks down. Nevertheless, we 
would expect the solutions to persist, with (j) (0) and ^ (tt) 
each moving toward one of the vacua. When )3 reaches 
the next critical value, {N + l)(N-|-4), two new solu¬ 
tions, with N+1 oscillations about (j)iop, will appear, but 
the previous ones will remain. Thus, foTN{N+l )<P< 
{N + l)(N-|-4), we should expect to find solutions with 
k = 0,1,2,... ,N oscillations. We have confirmed these 
expectations by numerically integrating the bounce equa¬ 
tions for various values of the parameters; the solutions 
for a typical potential with j3 = 70.03 are shown in Eig.|3 

The fact that the number of solutions should increase 
with j3 is physically quite reasonable. One would expect 






























the minimum distance needed for an oscillation about 
0 top, like the thickness of the bubble wall itself, to be 
roughly Hence, the number of oscillations that 

can fit on a sphere of radius H ^ should be of order 
^ Vl^"l In particular, this suggests that for 

j3 < 4 there should not even be a k = 1 Coleman-De 
Luccia bounce EH. 

It is thus somewhat puzzling to note the implications 
of Eq. for the case where A is negative. Here, in¬ 
creasing j3 through a critical value causes two solutions 
to merge into the Hawking-Moss solution and disappear, 
suggesting that the number of solutions is a decreasing 
function of j3. The resolution to this can be found by an¬ 
alytically and numerically examining various potentials 
that are unusually flat at the top. In all the cases we have 
examined, the number of solutions is governed by a pa¬ 
rameter 7 that measures an averaged value of \V"\/H^ 
over the width of the potential barrier. When 7 is suffi¬ 
ciently small, there are no bounce solutions (other than 
the Hawking-Moss, which is always present). As 7 is in¬ 
creased, new solutions appear at critical values. These 
first appear as solutions with finite values of 0 ( 0 ) — 0top- 
They then bifurcate, with 0 (0) for one solution moving 
toward a vacuum and 0 ( 0 ) for the other moving toward 
0top, eventually reaching it and disappearing when j3 is at 
a critical value. The net effect is that the number of solu¬ 
tions generally increases with 7 , although it is not strictly 
monotonic. 


INTERPRETING THE OSCILLATING 
BOUNCES 

How should these oscillating bounce solutions be in¬ 
terpreted? For the flat space bubble, a spacelike slice 
through the center of the bounce gives the initial con¬ 
ditions for the real-time evolution of the system after nu- 
cleation; these predict a bubble wall with a well-defined 
trajectory and a speed that soon approaches the speed 
of light. The interpretation of the Coleman-De Luccia 
bounce is similar. The main new feature here is the fact 
that the spacelike slice is finite. Formally, this corre¬ 
sponds to the fact that de Sitter spacetime is a closed 
universe, even though we expect the bubble nucleation 
process to proceed similarly in a spacetime that only ap¬ 
proximates de Sitter locally. 

The Hawking-Moss solution can be interpreted as cor¬ 
responding to a thermal fluctuation of all of de Sitter 
space (or, more plausibly, of an entire horizon volume) 
to the top of the potential barrier. Strictly speaking, clas¬ 
sical Lorentzian evolution would leave 0 at the top of the 
barrier forever. However, this is an unstable configura¬ 
tion, and so would be expected to break up, in a stochas¬ 
tic fashion, into regions that evolve toward one vacuum 


or the other. 

The oscillating bounce solutions yield a hybrid of 
these two extremes. The endcap regions near ^ —Q and 
^ = ^max clearly evolve into vacuum regions analogous 
to those from the Coleman-De Luccia bounce, while the 
intermediate, “oscillating”, region is like that emerging 
from a Hawking-Moss mediated transition. As with the 
Hawking-Moss solution, the bounce carries no informa¬ 
tion about the initial state, and there is not even any cor¬ 
relation between the vacua in the endcaps and the initial 
vacuum state. Thus, like Hawking-Moss, it is reminiscent 
of finite temperature tunneling in the absence of gravity, 
and provides evidence of the thermal nature of de Sitter 
spacetime. 

The relative importance of the various solutions de¬ 
pends on the values of their Euclidean actions. Although 
the details vary with the particular form of the poten¬ 
tial, the various regimes are characterized by a parame¬ 
ter 7 measuring an averaged value of |y"|///^. If 7 ^ 1, 
there is a Coleman-De Luccia bounce, a Hawking-Moss 
solution, and many oscillating bounces. However, the 
Coleman-De Luccia bounce has a much smaller action 
than the others, and so dominates. This is a regime of 
quantum tunneling transitions followed by determinis¬ 
tic classical evolution. At the other extreme is the case 
where 7^1, where the Hawking-Moss is the only solu¬ 
tion to the bounce equations. This is a regime of thermal 
transitions followed by stochastic real-time evolution. In 
between is a transitional region, with thermal effects still 
important. It is here that the oscillating bounces are most 
likely to play a role. 
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